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We evolve the network of global strings in the matter-dominated universe 

ON ' by means of numerical simulations. The existence of the scaling solution is 

^J^ ' confirmed as in the radiation-dominated universe but the scaling parameter 

£ takes a slightly smaller value, £ ~ 0.6 ±0.1, which is defined as £ = p s t 2 / ' \i 

Qh. with p s the energy density of global strings and \x the string tension per 

r- 1 \ unit length. The change of £ from the radiation to the matter-dominated 

universe is consistent with that obtained by Albrecht and Turok by use 

of the one-scale model. We also study the loop distribution function and 

find that it can be well fitted with that predicted by the one-scale model, 

where the number density n\{t) of the loop with the length I is given by 

m(t) = v/[t 2 {l + Kt) 2 } with v ~ 0.040 and « ~ 0.48. Thus, the evolution 

of the global string network in the matter-dominated universe can be well 

described by the one-scale model as in the radiation-dominated universe. 
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Cosmic strings are formed in a class of cosmological phase transitions [|IJ . Among them, 
gauged or local strings have been extensively studied because they could produce density 
fluctuations responsible for the large-scale structure formation and the cosmic microwave 
background anisotropy 0. A lot of analytical and numerical studies have been done to 
confirm that after some relaxation period a local string network enters the scaling regime. 
In this regime, the large-scale behavior of a network scales with the Hubble radius and 
the energy density of a network is given by 

Artring = i^/t 2 , (1) 

where £ is a constant irrespective of cosmic time. @-§|. Long strings intercommute to make 
loops [0-01, which decay through radiating gravitational waves |J. But an alternative 
scenario was recently proposed that long strings directly emit massive particles and lose 
their own energy |J. Thus, though the mechanism to drive a network to the scaling 
regime is still in dispute, the existence of the scaling solution is undoubted. 

On the other hand, global strings have been less investigated and considered only in 



the context of axion cosmology [1C 15]. It was assumed without direct verification that a 
global string network relaxes into the scaling solution just as the local counterpart, where 
long strings intercommute to create loops with a typical size that is comparable with the 
horizon scale and loops decay through emitting Nambu-Goldstone (NG) bosons. However, 
since global strings have a long-range interaction their dynamics could be different from 
that of local strings. In fact our numerical analysis of the evolution of a complex scalar 
field in 2+1 dimensions revealed that the global "strings" do not relax into the scaling 
solution but the number of defects per horizon volume increases logarithmically with 
time due to the long-range interaction [[U|. In the previous paper |L7|| we examined the 



evolution of a global string network in the radiation-dominated universe by numerically 
solving the equation of motion for a complex scalar field representing a global string. 
Then the above picture was confirmed. The scaling parameter £ takes a constant value 
and becomes £ ~ (0.9 — 1.3), irrespective of the cosmic time. Furthermore, one of the 
authors (M.Y.) |ISf has found that the number density rii(t) of the loop with the length 
/ can be well fitted by the formula predicted from the so-called one-scale model, that is, 
m(t) = z//[£ 3 / 2 (/ + Ktf/ 2 } with v ~ 0.0865 and k ~ 0.535. Thus, the evolution of a global 
string network in the radiation dominated universe is well described by the one-scale 
model proposed by Kibble [|T9"|-|2"T[ . 

In this paper, for a complete understanding, we study the evolution of a global string 
network in the matter-dominated universe, concretely, whether it can be well described 
by the one-scale model as in the radiation-dominated universe and, if at all, we clarify 
the relation between scaling parameters in both eras. 

We consider the following Lagrangian density for a complex scalar field $(x), 

£[$] = g„ v (d^)(d»$y - V eS [$,T] , (2) 

where g^ v is identified with the Robertson-Walker metric and the effective potential 
V eff [$,T] is given by 

V ee [$, T] = ^A($$t - v 2 ) 2 + l -\T 2 $$t , (3) 

which represents a typical second-order phase transition and the U(l) symmetry is broken 
below the critical temperature T c = \/3?7- 



For cosmological purposes, it would be desirable to trace the evolution of strings in the 
transition regime from the radiation-dominated era to the matter-dominated era. Due 
to the limitation of our computer powers, however, we concentrate on the evolution of a 
global string network during the matter-domination alone in this paper. Since the scaling 
property is expected to be reached irrespective of initial conditions if at all, we start 
simulations from a symmetric state with the equations of motion given in the matter 
domination. Strings are formed soon and they evolve in the matter-domination. As will 
be shown later, we confirm the scaling behavior in the matter-dominated regime and find 
that £ in this era is not different from that in the radiation-domination by more than 
a factor of 2. Hence we expect that transition from the radiation-dominated era to the 
matter-dominated era does not give rise to any significant cosmological effects, justifying 
our approach. 

In the matter-dominated universe, the equation of motion is given by 

*(x) + ZH*(x) - ^pV 2 $(x) = -KsM , (4) 

where the prime represents the derivative d/d& and R(t) is the scale factor which grows 
in proportion to £ 2 / 3 . We define a(T) [a(T) > 1] as a(T) = p ma ,t(T) / p Ta d(T) with p ma t(T) 
being the contribution to the energy density from non-relativistic particles and p ra d(T) 
being the contribution from relativistic particles at the temperature T. Then, 

a(T) = a c (T c /T) , (5) 

with a c = Pmak(T c ) / Piad(T c ) . Therefore, the Hubble parameter H = R(t)/R(t) and the 
cosmic time t are given by 

r, . . 87T 71 A 

H=a ^T^ T (6 > 
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where m p \ = 1.2 x 10 19 GeV is the Plank mass and g* is the total number of degrees of 
freedom for the relativistic particles. We define the dimensionless parameter ( as 

C ^7^ = [VSaJg^J ■ (9) 

In our simulation, we take ( = 8 and 4 to investigate ( dependence on the result. We take 
the initial time ti = t c /(2y / 2) corresponding to Ti = 2T C and the final time tf = 70(140) ti, 
where t c is the epoch with T = T c . Since the U(l) symmetry is restored at the initial 
time t = ti, we adopt as the initial condition the thermal equilibrium state with the mass 
squared, 



2 d 2 v cS m,T] 

m = 77--^ 



(10) 
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which is the inverse curvature squared of the potential at the origin at t = tj. 

Below we measure all of the physical quantities in units of tj. Then the equation of 
motion is given by 

*(») + f*(») - ^v»(.) = - (w + ^ - ^g=) •• , (id 

where A is set to unity for brevity. The scale factor R(t) is normalized as R(l) = 1. 

Using the second-order leap-frog method (see Ref. |T^] for details), we evolve the global 
string networks in the matter-dominated universe. In order to judge whether the global 
string network relaxes into the scaling regime in the matter-dominated universe, we give 
time development of the scaling parameter £ defined as p = C,p/t 2 . In our simulations, 
a lattice is identified with a part of a string core if the potential energy density there 
is larger than that corresponding to the field value of a static cylindrically-symmetric 
solution at r = 5x p h ys /\/2- 

We perform the simulations in seven different sets of lattice sizes, spacings, and ( (see 
Table H, |H]). In cases (1) and (5), the box size is nearly equal to the horizon volume 
(H^ 1 ) 3 and the lattice spacing to a typical width § ~ 1.0/(y/2r)) of a string at the final 
time tf. For each case, we simulate the system from 10 (Eqs. (l)-(6)) or 300 (Eq. (7)) 
different thermal initial conditions. 

If the simulation box is much larger than the horizon volume, it is reasonable to 
think that the boundary effect is negligible. But, in our simulation, the simulation box 
is comparable with or at most 4 3 times as large as the horizon volume at the final time 
of the simulation so that we should be careful to avoid possible boundary effects. In 
fact, a long-range force works between global strings so that the boundary effect cannot 
necessarily be neglected if the number of long strings in the simulation box is very small. 
As shown later, there are only a few long strings in our case so that the boundary effect 
could be significant. Hence in this paper we run simulations with two different boundary 
conditions of distinct features and employ a large enough simulation box so that the 
results with the different boundary conditions converge to each other. We thereby obtain 
a result free from the boundary effects. 

First, we adopt the periodic boundary condition, under which a string feels an at- 
tractive force from the boundary and there exists no infinite string so that strings can 
completely disappear in the simulation box if the Hubble radius becomes larger than 
the dimension of the simulation box. Figure [I] represents time development of £ with 
C = 8 (cases(l),(3),(4)) under this boundary condition. In case (1) with the smallest box, 
the boundary effect is so significant that strings tend to disappear, which is inconsistent 
with the result in Refs. (T^-^TJ . The larger the box size is, the less important the bound- 
ary effect is. In case (3) corresponding to the largest box simulations, the boundary effect 
is less significant so that £ seems to relax to a constant with £ ~ 0.49 ± 0.02. 

Next we adopt the reflective boundary condition, where V 2 $(x) on the boundary 
points disappears. Under this boundary condition, a string suffers a repulsive force from 
the boundary so that a string near the boundary intercommutes less often than that 
near the center of the simulation box because the partner to intercommute only lies in 
the inner direction of the boundary. Thus, the number of the strings tends to be more 
than that in the real universe. Figures |2| and |3] represent time development of £ with 
( = 8 [cases(l)-(4)] and ( = 4 [cases(5)-(7)], where £ becomes a constant irrespective of 



TABLE I. Three different sets of the simulations under the periodic boundary condition. 



Case 


Lattice 
number 


Lattice spacing 
[unit = tiR(t)\ 


c 


Realization 


Box size/H 1 
(at final time) 


e 


(1) 
(3) 
(4) 


128 3 
256 3 
256 3 


v/70/85 

v/70/85 

2^70/85 


8 
8 
8 


10 
10 
10 


l(at 70) 
2 (at 70) 
4(at 70) 


Disappearance 
0.50 ± 0.02 
0.49 ± 0.02 



TABLE II. Seven different sets of the simulations under the reflective boundary condition. 



Case 


Lattice 
number 


Lattice spacing 
[unit = UR(t)} 


c 


Realization 


Box size/H 1 
(at final time) 


£ 


(1) 


128 3 


v 7 ^/^ 


8 


10 


l(at 70) 


1.93 ±0.04 


(2) 


128 3 


2^70/85 


8 


10 


2 (at 70) 


1.03 ±0.03 


(3) 


256 3 


v 7 ^/^ 


8 


10 


2 (at 70) 


1.13 ±0.02 


(4) 


256 3 


2^/85 


8 


10 


4(at 70) 


0.72 ± 0.03 


(5) 


128 3 


^140/85 


4 


10 


l(at 140) 


1.36 ±0.06 


(6) 


256 3 


^140/85 


4 


10 


2(at 140) 


1.21 ±0.02 


(7) 


128 3 


2^70/85 


4 


300 


2 (at 70) 


1.25 ±0.02 
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FIG. 1. Time development of £ under the periodic boundary condition. Symbols (□) repre- 
sent time development of £. The vertical lines denote a standard deviation. 
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FIG. 2. Time development of £ in the cases from (1) to (4) with £ = 8 under the reflective 
boundary condition. Symbols (□) represent time development of £. The vertical lines denote a 
standard deviation over different initial conditions. 

time with (1) 1.93 ±0.04, (2) 1.03 ±0.03, (3) 1.13 ±0.02, (4) 0.72 ±0.03, (5) 1.36 ±0.06, 
(6) 1.21 ± 0.04, and (7) 0.72 ± 0.03. £ has larger values in smaller-box simulations due to 
the boundary effect as explained above. Also, for larger-box simulations, £ depends very 
little on £. 

From the results of the largest-box simulations containing the largest number of the 
Hubble volume at the final epoch, we conclude that £ converges to a constant £ ~ 0.6 ±0.1 
irrespective of the boundary conditions. This result will be supported later by comparing 
it with our previous results in the radiation-dominated regime p7|J18| through an analytic 
model of Albrecht and Turok ||. 

Next we investigate the loop distribution, which is predicted by Kibble's one-scale 
model fll9"pf2l|l as 



m(t) = v/[t 2 (l ± Ktf 



(12) 



where v is a constant, I is the length of a loop, and the logarithmic dependence of jjl is 
neglected. In contrast with local strings, the dominant energy loss mechanism of global 
strings is the radiation of the associated Nambu-Goldstone field [pTO]] . We define the 
radiation power P as P = Kfi where k is a constant. 

We determine whether the loop distribution in the simulation coincides with the above 
function. Since case (4) with the largest box takes too much time for one realization, we 
investigate the loop distribution for case (7) under the reflective boundary condition.[] 



1 £ t in case (7) is about twice as large as that in case (4) corresponding to the largest box 
simulation so that v in the real world may become half as large as that obtained in our simulation. 
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FIG. 3. Time development of £ in the cases from (5) to (7) with £ 
boundary condition. 



4 under the reflective 



The loop distribution is depicted in Fig. f| in case (7) at t = 40, 50, 60, and 70. Since long 
strings are rare, we cut the length of loops into bins with the width 5x5x. Then, we divide 
300 realizations into 6 groups comprised of 50 realizations and we summed the number 
of loops over 50 realizations for each groups. The dot represents the number of loops 
averaged over six groups and the dashed line represents the standard deviation. They 
can be simultaneously fitted with the above formula if one takes v ~ 0.040 and k ~ 0.48. 
Fittings for k and v are also given in Fig. [5| Thus, the loop production function as well 
as the large scale behavior of the string scales together for the global string network. Note 
that k takes almost the same value both in the matter- and radiation-dominated universe, 
which implies that NG bosons emission is not affected by the background universe and 
the large-scale behavior of a string network but is decided by the fundamental physics 
near a string segment. 

As discussed in 0, the one-scale model can predict the scaling parameter £ m in the 
matter domination from £ r in the radiation domination. From the Nambu-Goto action 
together with the intercommutation effect, the evolution of the string energy density p in 
the expanding universe is given by 



dp/dt = -3Hp + (l-2V 2 )Hp- cp/L, 
= -3Hp + H 2 Lp - cp/L, 



(13) 



where L = t/ \f^ is the characteristic scale with p = p/L 2 , V 2 is the average velocity 
squared of a string, and c is a constant representing the intercommutation efficiency.^ 



2 Though c may be dependent on V, we set c to be a constant as the zeroth order approximation. 







loop distributions 






8 


' A 








i i i i | 


; 


6 

4 


A 










: _ 


o 2 
p 10 


: t = 4o 




'. t = 50 " 














\ 


6 


i\ 




1 V ' ~~ 


simulation 
- fitting 


: 


2 


; t = 6o 




; t = 70 


'■| ,'"i , I 1 -,: 


■ I 



20 40 60 



100 20 40 60 80 100 



length \l/5x . 

FIG. 4. Loop distributions at t = 40, 50, 60, 70 are depicted. The number is summed over 
the box size [128((5x) 3 ] and 50 realizations for each groups. Bins are cut every 5x5x. 




FIG. 5. Fittings for k and v. The cross point represents the best fit values for k and v. The 
solid circle denotes 68% C.L. 



Defining 7 as 7 = H 1 /L, the evolution of 7 is given by 

^ = -|{c 7 2 -[2(ij"-i)-3] 7 -l}. (14) 

Then, in the radiation-domination with {H -1 ) = 2, we obtain the fixed point 7 r given by 

7r = (1 + yT+Ac)/2c. (15) 

On the other hand, in the matter domination with {H^ 1 ) = 3/2, we obtain the fixed point 
7 m given by 

7m = 1/vfc (16) 

Considering £ r = 7^/4 and £ m = 47^/9, we find the relation between £ r and £ m , 

e =- gr (it) 

U 9(2^+1)' U0 

Putting £ r ~ (0.9 — 1.3) |L8| into the relation, £ m is predicted to be 0.6 — 0.7, which 
coincides with our numerical results. 

In this paper, we investigated the evolution of a global string network in the matter- 
dominated universe. The network relaxes into the scaling regime as in the radiation- 
dominated universe but the scaling parameter £ m takes a smaller value, £ m ~ 0.6 ±0.1, 
which is consistent with the value predicted from £ r by use of the formula obtained by 
Albrecht and Turok |J. The loop distribution is also obtained and compared with that 
predicted by the one-scale model, where the number density ni(t) of the loop with the 
length I is given by ni(t) = v/{t 2 [l + nt) 2 }. With v ~ 0.040 and k ~ 0.48, the loop 
distribution function can be well fitted with that predicted by the one-scale model. Thus, 
the evolution of the global string network in the matter-dominated universe can be well 
described by the one-scale model as in the radiation-dominated universe. 

This work was partially supported by the Japanese Grant-in-Aid for Scientific Re- 
search from the Monbusho, Nos. 10-04558 (M.Y.), 11740146 (J.Y.), and "Priority Area: 
Supersymmetry and Unified Theory of Elementary Particles (# 707)" (J.Y. and M.K.). 
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